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Abstract 

A new (scalar) spectral decomposition is found for the Dirac system in two dimensions 
associated to the focusing Davey-Stewartson II (DSII) equation. Discrete spectrum in the 
spectral problem corresponds to eigenvalues embedded into a two-dimensional essential spec- 
trum. We show that these embedded eigenvalues are structurally unstable under small vari- 
ations of the initial data. This instability leads to the decay of localized initial data into 
continuous wave packets prescribed by the nonlinear dynamics of the DSII equation. 
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1 Introduction 



Gravity-capillary surface wave packets are described by the Davey-Stewartson (DS) system M 
which is integrable by inverse scattering tranform in the limit of shallow water || . In this paper, 
we study the focusing DSII equation which can be written in a complex form, 

iu t + u zz + u S g + A{g + g)u = 0, 

2g- z -(\u\ 2 ) z = 0, (1.1) 

where z — x + iy, z = x — iy, u(z, z, t) and g(z, z, t) are complex functions. This equation appears 
as the compatibility condition for the two-dimensional Dirac system, 

LfXz = -ULp 2 , Cp 2z = U(ft, (1.2) 

coupled to the equations for the time evolution of the eigenfunctions, 

iipit + ip lzz + mp 2z ~ u z ip 2 + Agipx = 0, 
- if2t + <f2zs + u z tpi - utp lz + Agip 2 = 0. (1.3) 

The DSII equation was solved formally through the d problem of complex analysis by Fokas 
and Ablowitz || and Beals and Coifman ||. Rigorous results on existence and uniqueness of 
solutions of the initial- value problem were established under a small-norm assumption p[. The 
small-norm assumption was used to eliminate homogeneous solutions of equations of the inverse 
scattering which correspond to bound states and radially symmetric localized waves (lumps) of the 
DSII equation. When the potential in the linear system becomes weakly localized (in L 2 but not 
in L 1 ), homogeneous solutions may exist and the analysis developed in Ref. |J is not applicable. 

The lump solutions were included formally in 0, where their weak decay rate was found, 
u ~ 0(i? _1 ) as R = \fx 2 + y' 2 — > oo. This result is only valid for complexified solutions of the 
DSII equation (when \u\ 2 is considered to be complex). The reality conditions were incorporated 
in the work of Arkadiev et al. where lumps were shown to decay like u ~ 0(i? -2 ). Multi-lump 
solutions were expressed as a ratio of two determinants B, or, in a special ratio of two 

polynomials but their dynamical role was left out of consideration. 

Recently, structural instability of a single lump of the DSII equation was reported by Gadyl'shin 
and Kiselev || §]. The authors used methods of perturbation theory based on completeness of 
squared eigenfunctions of the Dirac system [TU, |TT|]. A similar conclusion was announced by Yurov 



who studied Darboux transformation of the Dirac system |L2|, [L3|. 

In this paper, we present an alternative solution of the problem of stability of multi-lump solu- 
tions of the DSII equation. The approach generalizes our recent work on spectral decomposition of 
a linear time-dependent Schrodinger equation with weakly localized (not in L 1 ) potentials [T^| . We 



find a new spectral decomposition in terms of single eigenfunctions of the Dirac system. Surprisely 
enough, the two-component Dirac system in two dimensions has a scalar spectral decomposition. 
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In contrast, we recall that the Dirac system in one dimension (the so-called AKNS system) has a 
well-known 2x2 matrix spectral decomposition |[L5|| . 

Using the scalar spectral decomposition, we associate the multi-lump potentials with eigen- 
values embedded into a two-dimensional essential spectrum of the Dirac system. Eigenvalues 
embedded into a one- dimensional essential spectrum occur for instance for the time-dependent 
Schrodinger problem flj], Q . They were found to be structurally unstable under a small variation 
of the potential. Depending on the sign of the variation, they either disappear or become resonant 



poles in the complex spectral plane which correspond to lump solutions of the KPI equation [|14 

For the Dirac system in two dimensions, the multi-lump potentials and embedded eigenvalues 
are more exotic. The discrete spectrum of the Dirac system is separated from the continuous 
spectrum contribution in the sense that the spectral data satisfy certain constraints near the 
embedded eigenvalues. These constraints are met for special solutions of the DSII equation such 
as lumps, but may not be satisfied for a generic combination of lumps and radiative waves. As a 
result, embedded eigenvalues of the Dirac system generally disappear under a local disturbance of 
the initial data. Physically, this implies that a localized initial data of the DSII equation decays 
into radiation except for the cases where the data reduce to special solutions such as lumps. 

The paper is organized as follows. Elements of inverse scattering for the Dirac system are 
reviewed in Section 2, where we find that the discrete spectrum of the Dirac system is prescribed 
by certain constraints on the spectral data. Spectral decomposition is described in Section 3 with 
the proof of orthogonality and completeness relations through a proper adjoint problem. The 
perturbation theory for lumps is developed in Section 4 where some of previous results || || are 
recovered. Section 5 contains concluding remarks. Appendix A provides a summary of formulas 
of the complex (^-analysis used in proofs of Section 3. 

2 Spectral Data and Inverse Scattering 



Here we review some results on the Dirac system (|1.2| ) and discard henceforth the time dependence 
of u, g and (p. The potential u(z, z) is assumed to be non-integrable {u ^ L 1 ) with the boundary 
conditions, u ~ 0(|z| -2 ) as \z\ — > oo. 

2.1 Essential Spectrum of the Dirac System 

We define the fundamental matrix solution of Eq. ( |1.2j ) in the form |J, 

<p = [n(z, z, k, k)e ik \ X (z, z, k, k)e- fh "} , (2.1) 
where k is a spectral parameter, fi(z, z, k, k) and x{ z i z, k, k) satisfy the system, 

fi u = -ufi 2 , H2z = -ikfi 2 + Wt, (2.2) 

Xiz = ikxi~uX2, X2z = uxi- (2.3) 
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It follows from Eqs. (|2.2| ) and (|2.3|) that /j, and x are related by the symmetry constraint, 

X(z,z,k,k) = ap.(z,z,k,k), CT = ( 1 ' 
We impose the boundary conditions for fj,(z, z, k, k), 

lim u(z, z, k, k) — ei = ( \ ) . (2.5) 

|fc|^oo V V / 



Solutions of Eq. Q2.2|) with boundary conditions fl2.5|) can be expressed through the Green's 
functions as Fredholm's inhomogeneous integral equations |3|, |3j], 

^(z,z,k,k) = 1 ~~r f f dz Adz ( U fi 2 )(z', z'), (2.6) 

27U JJ z' — z 

» 2 (z,z,k,k) = f[%^(m)V,*)e-**'-'' ) ^* s -* ) . (2.7) 



Values of for which the homogeneous system associated to Eqs. (|2.6| ) and (|2.7| ) has bounded 
solutions are called eigenvalues of the discrete spectrum of the Dirac system. Let us suppose that 
the homogeneous solutions (eigenvalues) are not supported by the potential u(z,z). We evaluate 
the departure from analyticity of \x in the k plane by calculating the derivative dfi/dk directly 
from the system (H)-(|I3) 1 as 



^- = b(k,k)N,(z,z,k,k). (2.1 



Here b(k, k) is the spectral data, 



b(k, k) = ^- JJ dzAdz (Ufa) (z, z)e i{kz+kz \ (2.9) 

and N^z, z, k, k) is a solution of Eq. ( |2.2j) which is linearly independent of /J,(z,z,k,k) and 
satisfies the boundary condition, 

, lim N,(z, z, k, k) e ^ + ^ = e 2 = ( ° ) . (2.10) 

This solution can be expressed through the Fredholm's inhomogeneous equations, 

1 r r d.z! A dz' 

N lfX (z,z,k,k) = —jj——( u N 2ll )( Z ',z'), (2.11) 

N 2 Jz, z, k, k) = e -^ kz+k ^ + — [f dz ' Ad f ' (uN Xll )(z', ^ e -iKz-z')-fk(z-z')_ (212) 

2ni J J z' — z 
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The following reduction formula @j connects N^(^, z, k, k) and /j,(z, z, k, k), 

N M (z, z, k, k) = x(z, z, k, k) e - i{kz+ ~ k ~ z) = (Tfi(z, z, k, k)e~ i{kz+k ~ z) . (2.13) 

If the potential u(z,z) has the boundary values u ~ 0(|z|~ 2 ) as \z\ — > oo, (u L 1 ), then, 
the integral kernel in Eq. ( [2.9|) is not absolutely integrable, while Eqs. Q2.6|) and ( |2.7|) are still 
well-defined. We specify complex integration in the z-plane of a non-absolutely integrable function 
f(z, z) according to the formula, 



dz A dzf(z, z) = lim // dzAdzf(z,z). (2.14) 

R^oo J J\z\<R 

The same formula is valid for integrating eigenfunctions of the Dirac system in the k plane as 
well. In Section 3, we use fl2.14|) when computing the inner products and completeness relations 
for the Dirac system and its adjoint. 

2.2 The Discrete Spectrum 



Suppose here that integral equations Q2.6Q and (|2.7| ) have homogeneous solutions at an eigenvalue 
k = kj. The discrete spectrum associated to multi-lump potentials was introduced in Refs. [[|. |J. 
Here we review their approach and give a new result (Proposition 2.1) which clarifies the role of 
discrete spectrum in the spectral problem (|2.2|). 

For the discrete spectrum associated to the multi-lump potentials, an isolated eigenvalue k = kj 
has double multiplicity with the corresponding two bound states &j(z,z) and &j(z,z) 0. The 
bound state &j(z, z) is a solution of the homogeneous equations, 

<mm) = -hi S^f {u ^^ (2 - 15) 

<M*, 2) = 5^//^^(«*y)(^2')e-*^-* ,) -* ( *- - ' ) , (2.16) 
with the boundary conditions as \z\ — > oo, 

^(z,z)^^. (2.17) 



Equivalently, this boundary condition can be written as renormalization conditions for Eqs. ( |2.15| ) 
and (|2~T6D , 

— ffdz Adz(u$ 2i )(z,z) = 1, (2.18) 
2iri J J 

/ / dzAdzim^&zyto'+W = 0. (2.19) 
2ni J J 
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The other (degenerate) bound state <&j(z, z) can be expressed in terms of &j(z, z) using Eq. (|2.13|) , 

Qfa, z) = ct$ 3 (z, z)e~ i ^ z+ ~ k ^\ (2.20) 

The behaviour of the eigenfunction fj,(z, z, k, k) near the eigenvalue k = kj becomes compli- 
cated due to the fact that the double eigenvalue is embedded into the two-dimensional essential 
spectrum of the Dirac system ( |2.2j ). We prove the following result. 

Proposition 2.1. For smooth data b(k, k) e C 1 at k ^ kj, the eigenfunction fJt(z, z, k, k) has 
a pole singularity at k — > kj only if 

b = — f fdzAdz (zm X j) {z, z) e^ z+h ^ = 0. (2.21) 



Proof. Suppose fi(z, z, k, k) has a pole singularity at k — kj. Then, it can be shown from Eq. 
2~2] ) that the meromorphic continuation of fi(z, z, k, k) is given by the limiting relation, 



lim 



/j,(z,z,k,k) 



k-k 



'j J 



(Z + Zj)$j{z,z) +Cj&j{z,z), 



(2.22) 



where Zj, Cj are some constants. Using Eqs. ( |2.8|) , ( p.9|) , and ( p.!3|) , we find the differential 
relation for b(k, k), 



db _ b(k, k) 
dk ~ 2vr 



dz A dz (ufi 2 ) (z, z) — 



2ni 



dz A dz {zufix) (z, z)e i{kzJr ~ kt) . 



In the limit k — > kj, this equation reduces with the help of Eqs. (|2.18|) and (|2.22| ) to the form, 

db b(k, k) b 



dk k kj k kj 



where bo is given in Eq. ( [2.21| ). The reduced equation exhibits the limiting behavior of b(k, k) as 



k -> k 



3i 



b(k,k) 



-bo j j- In \k — kj\. 

rv rv>j 



(2.23) 



On the other hand, it follows from Eqs. fl2.13|) , ( p.20|) , and ( |2.22|) that N^(2;, z, k, k) has the 
limiting behavior, 

N M (z, z, k, k) - - Jl - (2.24) 



k kj 



According to Eqs. ( p. 23 ) and ( f2.24j ), the right-hand-side of Eq. ( |2.8| ) is of order O(bo\k—kj\ 1 ln|/c— 
kj\) as k — > kj. On the other hand, the left-hand-side of Eq. (|2.8| ) must be of order 0(1) in the 
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limit k kj according to Eq. ( |2.22|) . Therefore, the eigenfunction /j,(z, z, k, k) has a pole at 



k = kj only if the constraint 6 = holds. □ 



The limiting relation ( |2.22| ) was introduced by Arkadiev et al |J . However, the authors did not 
notice that the discrete spectrum is supported only by potentials which satisfy the additional con- 
straint (|2.21|) . In particular, such potentials include the multi-lump solutions for which b(k, k) = 
everywhere in the fc-plane. 



2.3 Expansion Formulas for Inverse Scattering 



Combining Eq. ( |2.8j ) for the essential spectrum and Eq. ( |2.22j ) for the discrete spectrum, we 
reconstruct the eigenfunction /j,(z, z, k, k) @, f|, 



fi{z,z,k,k) = ei + 2^ 



^ k kj QtTrt J J k k 



dk' A dk 1 

■b(k',k')N fl (z,z,k',k') 



(2.25) 



where n is number of distinct eigenvalues kj of double multiplicity. At k — > kj, this system is 
coupled with the algebraic system for the bound states, 



[z + Zj)&j(z, z) + CjQfo, z)=e 1 + J2 "t^u + i 



kn — ki 2ni J J k — k 



l+j ^3 ^ 



dk A dk 



b(k,k)N,,{z,z,k,k). (2.26) 



Expansion ( |2.25|) can be related to the inverse scattering transform for the potential u(z, z) 
R It follows from Eq. ( |2.2p that the eigenfunction n(z, z, k, k) has the asymptotic expansion 
as | A; | — > oo, 

/z^, z, fc, fc) = e a + ^^(z, z) + 0(|A;|- 2 ), (2.27) 



ik 



where /i2oo(z, z) = u(z,z) and 



Hioo(z, z) — — - — r 
Lixi 



dz' A cte' 



z — z 



(\u\*)(z',z f ). 



We deduce from Eqs. ( |2.25|) and ( p. 27] ) that the potential u(z,z) is expressed through the eigen- 
functions of the Dirac system in the form pf, 



u{z,z) = 



z,z) — 



3=1 



2tt 



dk A dk b(k, k) N 2fM (z, z, k, k). 



(2.28) 



Formulas ( ^.6[ ) to fl2.28| ) constitute a standard framework for the inverse scattering transform 
of the DSII equation with a new relation (|2 . 2 1 ) . The existence and uniqueness of solutions of the 
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Fredholm integral equations ( |2.6|) and (2.7) and the d problem Q2.8|) and ( |2.25|) were proved in |Q 
and || under the small- norm assumption for the potential u(x,y), 



sup \u\(x,y) ) [ I I \u(x, y)\dxdy ) < — . 



In this case n = and b(k, k) ^ 0. The nonlinear two-dimensional Fourier transform associated 
to this scheme was discussed in Examples 8-10 of Chapter 7.7 of Ref. ||17|| . Indeed, the connection 
formula (|2.28|) implies that there is a scalar spectral decomposition of u(z, z) through N 2fl (z, z, k, k) 



for n = 0. In order to close the decomposition, one could use Eqs. ( |2.9| ) and ( |2.13| ) to construct a 
"completeness relation" for the expansion of 5(z' — z) in the form, 



8(z' -z) = -—. JjdkA dkN 2 ^(z', z', k, k)N 2 ^(z, z, k, k). 

However, we show in Proposition 3.4 below that a completeness theorem for Eq. ( |2.2|) is different 
and is based on the set of eigenf unctions of the adjoint Dirac system. 



3 Basis for a Scalar Spectral Decomposition 



In this section, we specify the adjoint problem for the Dirac system ( |2.2|) and establish orthogo- 
nality and completeness relations. 

3.1 The Adjoint System 

The adjoint system for Eq. ( |2.2| ) is 

ML — ik\i\ — u^2, ttz = W$, (3- 1 ) 
which provides the balance equation, 

i(k' - k)rtfflto(k) = ^ ItfWtoik)] - ^ \tifflm(k)} . (3.2) 

The system ( |3.1| ) admits plane solutions fi a (z, z, k, k) and oscillatory- type solutions N^(z, z, k, k) 
with the boundary conditions, 

lim fi a (z,z,k,k) = e 2 , (3.3) 
lim N£(z, z, k, k)e- i{kz+ ~ k ~ z) = e v (3.4) 
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The adjoint eigenfunctions N^(z, z, k, k) can be expressed through the Green functions, 

N^{z,z) = e^ kz+h ^ - -L // dZ [ l ^\ uN^){z\z^ z - z '^ z -' z '\ (3.5) 

They are related to the adjoint eigenfunctions fi a (z, z, k, k) by the formula, 

N£(*, z, k, k) = -aji a {z, z, k, k)e i(kz+ ~ kz) . (3.7) 
Using this representation, we prove the following result. 

Lemma 3.1. The spectral data b(k, k) is expressed in terms of the adjoint eigenfunctions as 

b(k, k) = — JJdz A dz(uN^)(z, z). (3.8) 

Proof. Multiplying Eq. ( |3.5|) by u/ii(fc), integrating over dz A dz and using Eq. (|2.7| ), we 
express b(k, k) defined in Eq. (|2.9|) in the form, 



b(k, k) = — JJdz Adz [ufi!(k)N^(k) - ufi 2 (k)N^(k)\ . (3.9) 

On the other hand, multiplying Eq. ( |2.6| ) by uN^{k), integrating over dz A dz, and using Eqs. 
Q and (p|), we get Eq. (p|). □ 



Suppose now that k = kj is an isolated double eigenvalue of Eq. ( |2.2|) with the bound states 
3*j(z, z) and &j(z, z) given by Eqs. ( [2.15|) - Q2.20| ). Suppose also that k = k® is an eigenvalue of 
the adjoint system ( |3~T| ) with the adjoint bound states ^(z, z) and <frf(z,z). 

Lemma 3.2. If kj is a double eigenvalue of the Dirac system (|2.2p, then kj is also a double 



eigenvalue of the adjoint system (|3TT 



Proof. We use Eq. fl3.2| ) with = 3>j(z, z) and ^ a = &j(z, z) at k = kj and k! = kj and 
integrate over dz A <iz with the help of Eq. (A. 3) of Appendix A. The contour contribution of the 
integral vanishes due to the boundary conditions ( |2.17| ) and ( |3.1lQ and the resulting expression is 



(k] - kj) JJdzAdz ($^.$ 2j )(^, z) = 0. 

The relation kj = kj follows from this formula if the integral is non-zero at kj = kj (which is 
proved below in Eq. ( |3.22| )). The other possibility is when k a - ^ kj but is orthogonal to $2j- 
We do not consider such a non-generic situation. The other bound state at kj = kj can be 
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defined using the symmetry relation (see Eq. (|3.15|) below). □ 

The adjoint bound state (2, z) solves the homogeneous equations, 

2m 11 v — v 



z' — z 



with the boundary condition as \z\ — > oo, 



<f><*(z,z) 



e 2 



and the normalization conditions, 



1 



- — / dz A dz(u^ j )(z, z) = 1, 
2m J J 

dz Ndz{u^{z,z)e^ zVk ^ = 0. 



(3.10) 
(3.11) 

(3.12) 

(3.13) 
(3.14) 



In addition, the bound state &f(z,z) is related to &(z,z) according to the symmetry formula, 



(3.15) 



Using Eqs. (|3.1|)- Q3"TT5| ), we see that the adjoint eigenfunction fi a (z, z,k,k) satisfies relations 
similar to those for /j,(z, z, k, k), 



<9/x a 
~dk 



-b(k, k)N a Jz,z,k,k) 



and 



^lim 



V a (z } z,k,k)+ jV ) 



k-k 



3 J 



(3.16) 



(3.17) 



The expansions for inverse scattering transform of the adjoint eigenfunctions can be found in 
the form, 



^, *, *,*) = «*-£ ^fi£> -±-jj %^Hk\ k-)Nl(z, z, V, V) 

j = l 3 



(3.18) 



and 



(z'+ Zj )^(z', z r )-c^f(z', zf) = e 2 -g l ^ Z J^ ~JJ ^ffifo k)N«(z, z, k, k). (3.19) 
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3.2 Orthogonality and Completeness Relations 



Using the Dirac system fl2.2|) and its adjoint system fl3.1|) , we prove the orthogonality and com- 
pleteness relations for the set of eigenfunctions S = [N 2 ^(k, k), {$2j}"=i] and its adjoint set 
S a — [N^(k, k), {®ij}j =1 )- 

Proposition 3.3. The eigenfunctions N 2fl (z, z, k, k) and & 2 j{z, z) are orthogonal to the eigen- 
functions N° ' (z, z, k, k) and $ij(z, z) as follows 



(N^(k')\N 2 ,(k)) z = -27t 2 t5(k'-k), 
(N? fl {k)\&2j) z = (&ij\N 2 ^{k)) z = 0, 

(^ 2j ) z = 2m5 jh 

where the inner product is defined as 



(3.20) 

(3.21) 
(3.22) 



fo(W(*)>, 



dz A dz g(z, z, k', k') f(z, z, k, k). 



Proof. Using Eqs. fl2.12p and (^.5[), we expand the inner product in Eq. ( |3.20| ) as 

(N^(k')\N 2ll (k)) z = I 

+ JJdzAdz (uN^) (z, k'y-W^'^hiz) - [[dzA dz (uN lfi ) (z, k)e l{kz+k ~ z) h{z) 
// dz A dz (uNl) (z, k')e^ k ' z+k '^ 



d -^f (uN^) (z', k)e^ ,+k ^ [h(z) - h(z')) 



where 
and 



dz A dz e ^'-k)zM{k'-k)z = _2n 2 i8{k' - k) 



h{z) = — 
2m 



dz' A dz' 



z' — z 



,i(k'-k)z'+i(k'-k)z' 



i(k'-k)z+i(k'-k)z 



(3.23) 



(3.24) 



i(k' - k) 

The integrals Jo an d I\(z) axe computed in Appendix A. Using these formulas, we find the inner 
product in Eq. fl3.20| ) in the form, 



(Nl(k')\N 2fl (k)) z = -2ir*i5(k> -k) + — 
where the residual term TZ(k, k') is expressed in the form, 



K(k,k') 



K(k,k') 



dz A dz 



uN^N^-uN^N^k) 
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with the help of Eqs. fl2T|) and ([OP . We show that K(k, k') = by multiplying Eq. by 
uNin(k), integrating over dz A dz and using Eq. ( |2.11| ). 

The zero inner products in Eqs. (|3.21 ) and ( |3.22| ) for j ^ I are obtained in a similar way with 
the help of the Fredholm's equations for eigenfunctions <&j, N M , and N^. In order to find 
the non-zero inner product ( |3.22|) for j = I we evaluate the following integral by using the same 
integral equations, 



J J dz Adz $y/i 2 (fc) 



i(kj — k) 
1 

i(k — kj) 



dz A dz 



dz A dz 



Using Eq. (|3.13|) , the right-hand-side identifies to ^ h .y Substituting Eq. (|2.25| ) in the left-hand- 
side and using the zero inner products ( |3.21| ) and ( |3.22| ), we find Eq. ( |3.22| ) for j = I. □ 

Proposition 3.4. The eigenfunctions N 2fl (z, z, k, k) and Q 2 j(z,z) are complete with respect 
to the adjoint eigenfunctions (z, z, k, k) and $ij(z, z) according to the identity, 



8(z' - z) 



2n 2 i 



1 A 



dkAdkN?(z',z',k,k)N 2 ^z,z,k,k) --T i $Uz , ,z')$ 2j {z,z). (3.25) 

7T • -i 



i=i 



Proof. Using the symmetry relations ( |2. 13|) and (|3.7|) , we express the integral in Eq. (|3.25| ) as 



(N^(z')\N 2fl (z)) k = [[ dk A^^(^^^fc,% 1 (^,^,A ; ,fc)e ^fc ( 2 '- ^ ) +^ ^'- ^ ). (3.26) 



We use Eqs. ( [2.25| ), (|2.26|) , Q3.18| ), and (|3.19|) and find the pole decomposition for the integrand 
inEq. 



j=i K K j 



cj^Jz 1 , z')$'(z, z) + cMiz 1 , z')^(z, z) 



j=i k kj 



3=1 



(k-k 



3j 



+ 



2m 



dk' A dk' 
k' -k 



fl a 2 (z')bN^(z) - M^)bN^(z') (k\ k'). 



We substitute ( |3.27| ) into Eq. ( |3.26| ) and reduce the integral to the form, 



(3.27) 



(N^(z')\N 2ti (z)) k = Jo - 2ttE [^-(^^(M) + c^%{z' ,z')^{zrz)] h(k 3 
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- J2 z'^ijiz, z) [2n(z' - z)I\{kj) + J 2 (%)] 

3=1 



dkAdkh(k) Ji a 2 (z')bN ni (z) - Mz)bN^(z') (k,k)e 



-ik(z' —z)—ik(z' —z) 



(3.28) 



where the integrals Jo and are given in Eqs. (|3.23| ) and (|3.24j ) respectively, with z and 

interchanged, while the integral I 2 {kj) is defined by 



h(kj) = lim 



dk A d/c 



^ik(z' — z)+ik(z' — z) 



J J\k-kj\>e (k — kjY 



(3.29) 



The integral I 2 {k) is found in Appendix A in the form, I 2 {kj) = —2tt(z' — z)Ii(kj), such that the 
third term in Eq. ( |3.28| ) vanishes. In order to express the second term in Eq. ( |3.28| ) we use Eqs. 
(|2~26| ), ( gig ), (ET2CD , and fTPf) and derive the relation, 



z,z 



ikj (z 1 —z)+ikj (z' —z) 



3=1 



E *')<M*, f) + c^(z', z')& 2j (z, z) 

3=1 



n 1 " r r dk A dk r 



i=i 



27ri J J k — kj 

3=1 



<S> 2j (z)bN^(z') + ^{z')bN 2 ,{z)\ (k, k). 

(3.30) 



Using this expression and Eqs. (|3.23| ), Q3.24Q , and (|3.29| ), we rewrite Eq. ( |3.28|) in the form 



(N^(z')\N 2 ,(z)) k = -2nH6(z' - z) - 2m £ ^(z' ,z')$ 2j ( 

3=1 



Z,Z 



i f f - d 

+ —JJ dkAdk dk 



tfCO + E 



j=i k 



^(z) - E 



i$ 2j (z) 



j=i k kj 



The last integral vanishes according to Eq. (A. 3) of Appendix A and the boundary conditions 
Q and fl3j). □ 



Our main result for the spectral decomposition associated to the Dirac system ( |2.2| ) follows 
from the above orthogonality and completeness relations. 



Proposition 3.5. An arbitrary scalar function f(z,z) satisfying the condition f(z,z) 
0(|^| -2 ) as \z\ — > oo can be decomposed through the set S = [N 2fl (z, z, k, k), {$ 2 j(-2, z) }?=]]. 
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Proof. The spectral decomposition is defined through the orthogonality relations fl3.20|) -( j3.22|) 



as 



f(z, z) = JJdkA dka{k, k)N 2 ^(z, z, k,k) + aft 



2j\Z, Z), 



where 



a(k, k) 



4vr 2 



(N^(k)\f) z , ai = — (^|/> a 



(3.31) 



(3.32) 



Provided the condition on f(z, z) is satisfied, we interchange integration with respect to dz A dz 
and dk A dk and use the completeness formula (|3.25|) . □ 



The spectral decomposition presented here is different from that of Kiselev [K| |TT[. In the 
latter approach, the function f(z, z) is spanned by squared eigenfunctions of the original problem 
(|1.2j ) defined according to oscillatory-type behaviour at infinity. In our approach, we transformed 
the system (|1.2j ) to the form ( p.2|) and defined the oscillatory-type eigenfunctions according to the 
single eigenfunctions N^z, z, k, k). We also notice that the (degenerate) bound states ^(z,z) 
are not relevant for the spectral decomposition, although they appear implicitly through the 
meromorphic contributions of the eigenfunctions N At (z, z, k, k) at k = kj (see Section 4). 



4 Perturbation Theory for a Single Lump 

We use the scalar spectral decomposition based on Eq. ( |3.31| ) and develop a perturbation theory 
for multi-lump solutions of the DSII equation. We present formulas in the case of a single lump 
(n = 1), the case of multi-lump potentials can be obtained by summing along the indices j, I 
occuring in the expressions below. 

The single-lump potential u(z,z) has the form B, 



U{Z, z 



\Z + Zi 2 + \Cj 



i(kjZ+kjz) 



(4.1) 



where complex parameters. The associated bound states follow from Eqs. (|2.26|) and 

(|3~T9| ) as 



<$>j(z,z) 



\z + zA 2 + \ca 



\z + zA 2 + \ca 



Cj e 



Z + Zj 
z4 

i(kj z+kjz) 



*-i & %{kj z-\-kj z) 



Z + Zj 



(4.2) 
(4.3) 



We first consider a general perturbation to the single lump subject to the localization condition, 



Au ~ 0(|z| 2 ) as \z\ 
term Au(z, z). 



oo. We then derive explicit formulas for a special form of the perturbation 
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4.1 General Perturbation of a Single Lump 

Suppose the potential is specified as u e = u(z, z) + eAu(z, z), where u(z, z) is given by Eq. (|4.1| ) 
and Au(z,z) is a perturbation term. Two bound states $j(z, z) and &'Az,z) are supported by a 
single-lump potential u(z,z) at a single point k = kj. The spectral decomposition given by Eq. 
( |3.31| ) provides a basis for expansion of z, k, k) at k = k, 



^2\ Z i Z i K ) 



dk A dk a(k, fc)A^ 2M (z, z, k, k) + a.,^'^, z) 



(4.4) 



where a(/c, /c) and <x,- are defined by Eq. Q3.32| ) and depend on the parameter k. The other 
component n\(z, z, k, R) can be expressed from Eq. (|2.2|) as 



jJl^iyZ^ hzj 



dk A dk a(k, k)Ni^(z, z, k, k) + ctj&^^z, z) + eAfii(z, z), (4.5) 



where the remainder term A/j,i(z,z) solves the equation, 

( A ^i) 5 = -Aw/4 
We write the solution of this equation in the form, 



dz' A dz' 
z' — z 



{AuvQ {z', z! 



(4.6) 



subject to the boundary condition as \z\ — > oo, 

Afii(z,z) -> A + Oiz- 1 ), 

where A is an arbitrary constant. Using the explicit representation ( |4.6| ), we transform Eq. (fX 
into the system of integral equations for a(k, k) and a>j, 



a(k, k) 



An 2 i{k — k) 



dk' A dk' K(k, k, k', k')a(k', k') + Kj(k, k)a d + R(k } k)A 



0(e 2 ), 



a ; 



2n(kj — k) 



J J dkAdk Pj(k, k)a(k, k) + K jt ai + RjA 



+ 0(6 2 ), 



(4.7) 
(4.8) 



where 



K{k,k,k',k') = (N^AOIN^'))^, Kj(k,k) = (^(k)\^) Au , 

Pj(k, k) = (^|N Ai (A;)>Au, K fl = (*?|*i)a«, 
and the scalar product for squared eigenfunction is defined as p], O] 



dz A dz 



h(z, z)ft(z, z, k, k)gi(z, z, k', k!) + h(z, z)f 2 (z, z, k, k)g 2 (z, z, k', k' 
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The non-homogeneous terms R(k, k) and Rj can be computed exactly as 

R{k, k) = J J dzAdz (uN^) (z, z) = 2nb(k, k), 



Rq 



dz A dz [u^ljj (z, z) 



-27U, 



where b(k, k) = if n 7^ 0. We solve the system of equations ( |4.7| ) and ( [4.8| ) asymptotically for 
k = kj + eAft and An ~ 0(1). The leading order behaviour of the integral kernels follows from 
the asymptotic representation ( |2.24| ) as k — > kj, 



I^l (^k j k^ k . k 



K 



(k kj)(k' kj 



Kj(k,k) - %P " 



k kj 



Pj (k,k) iP " 



k kj 



(4.9) 



where 



Kjj = (<&f |<^>A„, Pjj = (*?\*j)Au, PjJ = -<*5I*}>A«. (4-10) 

Here we have used the symmetry constraints ( |2.20|) and ( p. 15 ). The leading order of a(k,k) as 
k — > kj follows from Eq. (|4.7| ) as 



a(k, k) 



2ir(k-K)(k- kj)' 

where j3j is not yet defined. We use Eq. (A. 4) of Appendix A to compute the integral term, 

dk A dk 2iri 



(k-K)(k-kjf 

and reduce the system of integral equations ( (4.7[) and 

— 27rAfta,- 



kj Kj 



to an algebraic system as k — > k 



j * 



2irARf3j 



~Pjj a j ~ KjjPj- 



(4.11) 
(4.12) 



If Pjj 7^ 0, the determinant of the above system is strictly positive. Therefore, homogeneous 
solutions at A = (bound states) are absent for e ^ 0. This result indicates that the double 
eigenvalue at k — kj disappears under a generic perturbation of the potential u(z, z) with Pjj ^ 
(see also Ref. |§). 

For A ^ 0, we find inhomogeneous solutions of Eqs. ( 4.11 ) and ( 4.12j ), 



2mA 


( K Kjj + 2vrA/^ 




\Kjj + 2-kAk 


2 + 


P 

r 33 


12 



ft 



-2vrzv4P, 



Kjj + 2ttA/s:| 2 + |P 7 -,-| 2 



(4.13) 



33 \ 
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The eigenfunction fj, e (z, z, k, k) given by Eqs. (4.4) and (|4.5|) satisfies the boundary condition 
if A = e _1 and has the following asymptotic representation, 



/X [Z, Z, K, Kj 



ei + 



2m 


2tt(k — kj) + eKjj 


<f>j{z,z) 


2meP^>(z,z) 


\2-k(k- kj) + eKjj\ 2 + \ePjj\ 2 


\M*-kj) + eK»\ 2 + \ePjj\ 2 



+ Afj, e {z,z) 



(4.14) 



where the term A[i e (z, z) is not singular in the limit e — > and K — > kj. 

In the limit e — > 0, K ^ kj, we find a meromorphic expansion for fi e (z, z, K, R) as 



/X (Z, 2, ft, Kj 



ei + n ! ; + e 



2m\K — kj) 2 



2ni\K — kj\ 2 



+ 0(6 2 



(4.15) 



It is clear that the double pole can be incorporated by shifting the eigenvalue kj to 



h e — h. - eKjj 
3 3 2tt 



The other double-pole term in the expansion ( |4.15| ) has a non-analytic behaviour in the /c-plane 
and leads to the appearance of the spectral data b e (K, R) = eAb(n, R) which measures the departure 
of ^ e (z, z, k, R) from analyticity according to Eq. (|2.8| ). We find from Eqs. (|2.9|) and ( |4.15|) that 
the spectral data Ab(n, R) has the following singular behaviour as k —>■ kj 



A6(k, k) 



P. 



■jj 



2tc\k — ki 



(4.16) 



Thus, if Pjj 7^ the analyticity of fj, e (z,z, k,k) is destroyed and the lump disappears. This 
conclusion as well as the analytical solution (|4.13|) agree with the results of Gadyl'shin and Kiselev 
|| H where the transformation of a single lump into decaying wave packet was also studied. 
In the other limit e ^ and k — > kj we find another expansion from Eq. (|4.14j) , 



/X [Z, Z, K, Kj 



ei 



2m 



eP ' J 

tr 33 



&j{z,z)+0{K-k*). 



We conclude that the eigenfunction fJ, e (z, z, k, k) is now free of pole singularities 
marize the main result in the form of a proposition. 



(4.17) 
I. We sum- 



Proposition 4.1. Suppose u(z,z) is given by Eq. Q4.1| ) and Au(z,z) satisfies the constraint, 

Pjj = <*?l*i>A« + 0. 



Then, the potential u e = u(z, z) + eAu(z, z) does not support embedded eigenvalues of the Dirac 
system (|2.2|) for e ^ 0. 
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4.2 Explicit Solution for a Particular Perturbation 

Here we specify Cj = ce l9 , where c and 9 are real, and consider a particular perturbation Au(z, z) 
to the lump u(z, z) ( |4.1|) in the form, 

Au{z,z) = Q(z,z)e i{k ^ z+ ^~ z+e) , 

where Q(z, z) is a real function. Using Eqs. ( |4.1(J| ) , and (|4.3|), we find explicitly the matrix 

elements Kjj and Pjj, 

^ — J J dz A ^ ^ffg^ - 1 // * A * + . 

The element Pjj can be seen as a correction to the field energy, 



N = - J J dz A dz\u e \ 2 (z, z) = N + iecPjj + 0(e 2 ) 



i 

2 

where Nq — ir is the energy of the single lump solution (independent of the lump parameters kj 
and Cj). Thus, a perturbation which leads to the destruction of a single lump, that is with Pjj ^ 0, 
changes necessarily the value for the lump energy iVo- 

5 Concluding Remarks 

The main result of our paper is the prediction of structural instability of multi-lump potentials 
in the Dirac system associated to the DSII equation. The multi-lump potentials correspond to 
eigenvalues embedded into a two-dimensional continuous spectrum with the spectral data b(k, k) 



satisfying the additional constraint (|2.21| ). In this case, there is no interaction between lumps 
and continuous radiation. However, a generic initial perturbation induces coupling between the 
lumps and radiation and, as a result of their interaction, the embedded eigenvalues disappear. 
This result indicates that the localized multi-lump solutions decay into continuous wave packets 
in the nonlinear dynamics of the DSII equation (see also Refs. || £5|). 

This scenario is different from the two types of bifurcations of embedded eigenvalues discussed 
in our previous paper ||14|| . The type I bifurcation arises from the edge of the essential spectrum 
when the limiting bounded (non-localized) eigenfunction is transformed into a localized bound 
state. The type II bifurcation occurs when an embedded eigenvalue splits off the essential spec- 
trum. Both situations persist in the spectral plane when the essential spectrum is one-dimensional 
and covers either a half-axis or the whole axis. However, in the case of the DSII equation, the es- 
sential spectrum is the whole spectral plane and embedded eigenvalue can not split off the essential 
spectrum. As a result, they disappear due to their structural instability. 
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Appendix A. Formulas of the ^-analysis 

Here we reproduce some formulas of the complex <9-analysis |L7 to compute the integrals Iq, Ii(z) 
and defined in Eqs. ( |3. 23| ) , ( |3.24| ), and ( |3.29| ). We define the complex integration in the 

z-plane by 

dz A dzf(z, z), 



dz A dzf(z, z 

where dz A dz = —2idxdy. The complex 5(z) distribution is defined by 

dz A dzf(z)5(z - zo) = -2if(z ), (A.l) 

where 5 (z) = S(x)5(y). In particular, the ^-distribution appears in the (^-analysis according to the 
relation [FTi 



d_ 

dz L z — z 

Computing the integral 1$, we get the formula, 



ii5(z - zo). 



(A.2) 



dz A dze 



— ikz+ikz 



-2% 



dxdye 



2iRe(k)x-2iIm(k)y 



-2w 2 i5(k), 



which proves the identity (gTgg ). 

Using the Green's theorem W%, one has the integration identity, 



(A.3) 



where D is a domain of the complex plane and C its boundary. The generalized Cauchy's formula 

f{z',z')dz> 



has the form 17 



or, equivalently, 



f(z,z) 



1 

27T? Jc Z' — Z 

1 



+ 



f{z',z')dz' + 1 



1 r r dz' A dz' df 
2ni J Jd z' — z dz' ' 

dz' A dz' df 



2ni Jc z' — z ' 2ni J Jd z' — z dz' 
In order to find the integral h{z) we use Eq. (A. 5) with 

f{ zrz ) = ^ kz+M \ MO 



(AA) 
(A5) 
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and choose the domain D to be a large ball of radius R (see Eq. (|2.14j) ). The boundary value 
integral vanishes since 



lim / ^l e i(kz+kz) = _ 27r - Um j^ 2 \k\R) = 0, {A.6) 

R-^oo J\z\=R Z R^oo 

where Jq(z) is the Bessel function. Equation (A. 5) for the function f{z,z) then reduces to Eq. 



In order to compute the integral ^(^o) > we a Pply Eq. (A. 3) with / x = and 

f 2 (z,z) = -J— e i ^+~^). 
z- z 



The domain D is chosen as above. The boundary value integral vanishes again, 

li m / <^ e i(kz+kz) = _ 2ni Um j_ 2 (2\k\R) = 0, (A.7) 

R->oo J\z\=R Z R^oo 

where J-2{z) is the Bessel function. Equation (A. 3) for the function f2(z, z) reduces to Eq. ( |3.29| ). 
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